Measurement of and v a induced neutral current single 7r° production cross sections 

on mineral oil at E v ~ 0(1 GeV) 

A. A. Aguilar-Arevalo, 11 C. E. Anderson, 19 A. 0. Bazarko, 16 S. J. Brice, 7 B. C. Brown, 7 L. Bugel, 5 J. Cao, 15 
L. Coney, 5 J. M. Conrad, 14 D. C. Cox, 10 A. Curioni, 19 Z. Djurcic, 5 D. A. Finley, 7 B. T. Fleming, 19 
R. Ford, 7 F. G. Garcia, 7 G. T. Garvey, 12 J. Gonzales, 12 J. Grange, 8 C. Green, 7 ' 12 J. A. Green, 10 ' 12 
T. L. Hart, 4 E. Hawker, 3 ' 12 R. Imlay, 13 R. A. Johnson, 3 G. Karagiorgi, 14 P. Kasper, 7 T. Katori, 10 ' 14 
T. Kobilarcik, 7 I. Kourbanis, 7 S. Koutsoliotas, 2 E. M. Laird, 16 S. K. Linden, 19 J. M. Link, 18 Y. Liu, 1 Y. Liu, 15 

W. C. Louis, 12 K. B. M. Malm, 5 W. Marsh, 7 C. Mauger, 12 V. T. McGary, 14 G. McGregor, 12 W. Metcalf, 13 
P. D. Meyers, 16 F. Mills, 7 G. B. Mills, 12 J. Monroe, 5 C. D. Moore, 7 J. Mousseau, 8 R. H. Nelson, 4 P. Nienaber, 17 
J. A. Nowak, 13 B. Osmanov, 8 S. Ouedraogo, 13 R. B. Patterson, 16 Z. Pavlovic, 12 D. Perevalov, 1 C. C. Polly, 7 
E. Prebys, 7 J. L. Raaf, 3 H. Ray, 8 ' 12 B. P. Roe, 15 A. D. Russell, 7 V. Sandberg, 12 R. Schirato, 12 D. Schmitz, 5 
M. H. Shaevitz, 5 F. C. Shoemaker*, 16 D. Smith, 6 M. Soderberg, 19 M. Soreb\ 5 P. Spentzouris, 7 J. Spitz, 19 
I. Stancu, 1 R. J. Stefanski, 7 M. Sung, 13 H. A. Tanaka, 16 R. Tayloe, 10 M. Tzanov, 4 R. G. Van de Water, 12 
M. O. WascW, 13 D. H. White, 12 M. J. Wilking, 4 H. J. Yang, 15 G. P. Zeller, 7 and E. D. Zimmerman 4 

(The MiniBooNE Collaboration) 

1 University of Alabama; Tuscaloosa, AL 354-87 
2 Bucknell University; Lewisburg, PA 17837 
s University of Cincinnati; Cincinnati, OH 45221 
4 University of Colorado; Boulder, CO 80309 
•' Columbia University; New York, NY 10027 
6 Embry Riddle Aeronautical University; Prescott, AZ 86301 
7 Fermi National Accelerator Laboratory; Batavia, IL 60510 
8 University of Florida; Gainesville, FL 32611 

9 University of Illinois; Urbana, IL 61801 
10 Indiana University; Bloommgton, IN 47405 
11 Instituto de Ciencias Nucleares, Universidad Nacional Autonoma de Mexico, D.F. 04510, Mexico 
12 Los Alamos National Laboratory; Los Alamos, NM 87545 
Louisiana State University; Baton Rouge, LA 70803 
' Massachusetts Institute of Technology; Cambridge, MA 02139 
15 University of Michigan; Ann Arbor, MI 48109 
Princeton University; Princeton, NJ 08544 
Saint Mary's University of Minnesota; Winona, MN 55987 
18 Virginia Polytechnic Institute & State University; Blacksburg, VA 24061 
19 Yale University; New Haven, CT 06520 

MiniBooNE reports the first absolute cross sections for neutral current single it production 
on CH2 induced by neutrino and antineutrino interactions measured from the largest sets of NC 
7T° events collected to date. The principal result consists of differential cross sections measured as 
functions of tt° momentum and ty° angle averaged over the neutrino flux at MiniBooNE. We find total 
cross sections of (4.76±0.05 s t a t±0.76 S i/s) X 10 -40 cm 2 /nucleon at a mean energy of {E„) = 808 MeV 
and (1.48±0.05 s tat ± Q.23 syB ) x 10 _ cm 2 /nucleon at a mean energy of {E v ) = 664 MeV for and 
induced production, respectively. In addition, we have included measurements of the neutrino 
and antineutrino total cross sections for incoherent exclusive NC l7r° production corrected for the 
effects of final state interactions to compare to prior results. 



I. INTRODUCTION 



Neutral current neutrino interactions producing a sin- 
gle 7T° (NC l7r°) constitute a substantial background for 
experiments searching for — > v e oscillations. NC l7r° 
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events are prone to mimicking single electrons — the sig- 
nature sought in such z/ e -appearance searches — because 
one of the two photons from the 7r° decay may escape 
detection. In MiniBooNE, NC l7r° production poses one 
of the largest backgrounds: it is second only to events 
induced by intrinsic v e in the beampQ. As such, abso- 
lute measurements of NC l7r° production at energies of 
0(1 GeV) are crucial to constraining this background, es- 
pecially as it applies to future long-baseline experiments. 

A measurement of NC l7r° production can also be used 
to test and refine models of single 7r° production, which 
vary widely in their predictions at these energies [2lll7|. 
These models categorize exclusive NC l7r° production on 
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nuclei by final state as either coherent or incoherent. Pro- 
duction leaving the nuclear target in the ground state 
is defined as coherent, otherwise it is defined as inco- 
herent. Prior measurements of NC 1tt° production were 
typically limited in scope, having addressed incoherent 
and coherent production separately, and suffered from 
low statistics. The earliest results were total cross sec- 
tions measured as ratios normalized to various charged 
current pion production channels [T%l - f2"2"] . Later, stud- 
ies of absolute NC 1tt° production were performed. Ab- 
solute measurements of incoherent NC lir° production 
were reported by Aachen-Padova[2"3"] (albeit in a foot- 
note) and in a more recent reanalysis of Gargamelle 
data [H], both at neutrino energies near 2 GeV. The dis- 
tinct signature of coherent NC lir° production — a for- 
ward emitted ir° and a target left in its ground state — 
permits absolute measurements of coherent NC lir° pro- 
duction. These measurements were carried out under 
a variety of circumstances [531 I2"5l428| . While measure- 
ments regarding such exclusive production are valuable, 
the total yield of NC l7T° production is often more im- 
portant to modern-day neutrino oscillation experiments. 
To address this need, inclusive NC lir° and NC w° mea- 
surements, reported as flux-averaged cross section ratios 
relative to CC production, have been recently performed 
by K2K[1H] and SciBooNE[3D], respectively. Collectively, 
prior experiments have recorded a few thousand neutrino 
and a few hundred antineutrino NC lir° interactions. 

In this paper, MiniBooNE reports the first measure- 
ments of absolute inclusive NC 1tt° cross sections (not 
normalized as ratios) for both neutrino and antineutrino 
scattering. We define signal NC lir° events to be NC in- 
teractions wherein only one ir° and no additional meson 
exits the target nucleus (no requirement on the number 
or identity of outgoing nucleons is made) . This definition 
is consistent with that used at K2K[29 . It is specifically 
chosen because final state interactions (FSI) dramatically 
alter the experimentally observed products of the origi- 
nal neutrino interaction on a nuclear target, but are not 
well understood. As particles in the final state transit 
the nucleus, they can scatter, be absorbed, or undergo 
charge exchange. The observation of NC l7r° interactions 
in an experiment will be depleted by the effects of absorp- 
tion and charge exchange (7r°p — > 7r + n, 7r°n — > 7r~p); 
however, it can also be enhanced by additional chan- 
nels entering the sample if a ir° is produced via FSI (e.g. 
7r + n — » 7r°p, 7r~p — > 7r°n, or 7T° production from nucleon 
rescattering). Ultimately, it is this observed rate of 7r° 
production, regardless of the initial interaction, that is 
relevant to neutrino oscillation experiments operating on 
nuclear targets. Hence, the definition of our signal, one 
constructed in terms of the observed final state, directly 
addresses the requirements for — > v e oscillation exper- 
iments. At the same time, the inclusivity of the defini- 
tion reduces the dependence of the measurement on the 
assumed models of FSI and single 7r° production. Here- 
after, we use "NC l7r°" to refer to this inclusive definition 
unless explicitly stated otherwise. Under this definition 



and in a calculated effort to reduce model dependence, 
we present the first absolute differential and total cross 
sections for and P p induced NC l7r° production; the 
interactions occurred on CH2. Since the neutrino en- 
ergy cannot be measured for each interaction, the cross 
sections are necessarily averaged over the neutrino flux 
at MiniBooNE. Specifically, we have measured cross sec- 
tions as a function of ir° momentum (p^o) and 7r° angle 
relative to the interacting neutrino (cos#„.o). Together, 
these measurements can yield important information on 
FSI effects, which are a strong function of 7r° momentum, 
and the production mechanism (coherent versus incoher- 
ent), which is a strong function of ir° angle. 



II. THE EXPERIMENT 

MiniBooNE receives neutrinos from the Booster Neu- 
trino Beam at Fermilab. 8 GeV protons extracted from 
the Booster synchrotron are delivered to a beryllium 
target; neutrinos result from the decays of secondary 
mesons produced by interactions in the target. The tar- 
get is housed in a magnetic horn which focuses charged 
mesons of a selected sign and defocuses mesons of the 
opposite sign. A beam which is predominately com- 
posed of either neutrinos or antineutrinos can be pro- 
duced by choosing the polarity of the horn current. In 
neutrino mode, with a mean energy of 808 MeV com- 
prise 93.6% of the flux and contamination from v^, v e , 
and v e comprise 5.86%, 0.52%, and 0.05% of the flux, re- 
spectively. Wrong-sign [50J (WS) contamination impacts 
the antineutrino mode flux to a greater degree. In an- 
tineutrino mode, with a mean energy of 664 MeV 
comprise 83.73% of the flux and contamination from v^, 
24, and P e comprise 15.71%, 0.2%, and 0.4% of the flux, 
respectively (31]. 

The detector[32] consists of a 12.2 m diameter spher- 
ical vessel filled with 818 tons of undoped mineral oil 
situated 541 m from the target. The containment vessel 
is segmented by an optical barrier into a 5.75 m radius 
inner tank region and an additional 0.35 m veto region. 
The surface of the inner tank is instrumented with 1280 8- 
inch photomultiplier tubes (PMTs), which provide 11.3% 
photocathode coverage. The tank PMTs capture the pat- 
tern of light generated by charged products of neutrino 
interactions. Particles above Cherenkov threshold emit 
directional light conically about the particle track which 
produces a ring on the tank surface. Isotropic scintilla- 
tion light emitted by certain constituents of the mineral 
oil is also detected by the PMTs. The veto region, which 
is instrumented with 240 PMTs, is used to detect light 
due to particles entering or exiting the detector. 

Neutrino interactions in MiniBooNE are simulated us- 
ing the v3 Nuance event generator[33] coupled to a 
GE ANT3-based [34] detector Monte Carlo. Single n° pro- 
duction is predicted according to the models of Rein and 
Sehgal (R-S)[H[S] as implemented in Nuance with two 
exceptions. First, we modify Nuance to incorporate 
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Multi-7r 
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17% 
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DIS 


< 1% 
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NC TT* 


2% 


2% 


K, p, »7 Prod. 


< 1% 


< i% 



TABLE I: Predicted fractional composition of NC Itt si K - 
nal events in neutrino and antineutrino modes broken down 
according to exclusive channel at the neutrino interaction ver- 
tex. 



nonisotropic A decays. Second, the relative contribution 
of coherent and incoherent exclusive NC l7r° production 
is further adjusted using a prior measurement [35] : coher- 
ent pion production is reduced by 35% and incoherent is 
increased a corresponding 5% to preserve total ir° pro- 
duction. The FSI model in Nuance accounts for the 
rescattering of all hadrons during nuclear transit; the 
pion absorption factor described in the R-S model of 
coherent pion production is omitted in lieu of the Nu- 
ance FSI model. In all, we predict 94% of observed 
NC l7r° production to involve the production of a ir° at 
the neutrino interaction vertex; the fraction rises to 97% 
in antineutrino mode. A breakdown of the composition 
of NC l7r° production by exclusive interaction channel is 
listed in Table [I] The R-S models predict a smaller in- 
coherent pion production cross section for antineutrinos 
than for neutrinos, but similar coherent pion production 
cross sections for both. As a result, the Monte Carlo 
predicts that the fraction of NC l7r° production that is 
coherent pion production is larger in antineutrino mode 
than in neutrino mode. In principle, this effect makes an- 
tineutrino scattering more sensitive to the coherent pion 
production mode. The mineral oil target, which consists 
largely of long alkanes and cycloalkanes, is simulated as 
CH 2 in Nuance. 21% of NC l7r° production is predicted 
to occur on free nucleons (hydrogen). This fraction is 
greater than the fraction of nucleons in CH 2 (14.3%) be- 
longing to H because nuclear effects (predominately pion 
absorption) diminish the cross section on carbon. 

III. SELECTION AND RECONSTRUCTION 

Before events are reconstructed, a series of simple cuts 
are made. Events are decomposed into sets of PMT hits 
clustered in time (subevents). Selected NC lir° candi- 
dates are required to have (1) only one subevent and 
that subevent is coincident with the 1.6/zs neutrino beam 
pulse. Multiple subevents arise principally from muon 
decays — a signature of charged current events or ir pro- 
duction. Further cuts require that the single subevent 
possess (2) fewer than 6 PMT hits in the veto region and 
(3) greater than 200 PMT hits in the tank region. The 
veto hits requirement removes uncontained events as well 
as events with particles entering the detector during the 
beam pulse. The tank hits requirement reduces the con- 



tamination from NC elastic events and eliminates events 
containing a decay electron from a cosmic muon entering 
the tank before the beam. 

After the preliminary cuts, the remaining events are 
reconstructed in order to measure kinematic variables 
and perform particle identification. The reconstruc- 
tion algorithm takes the form of a track-based, least 
negative-log-likelihood fit performed under various par- 
ticle hypotheses [3"rj]. Four hypotheses are used in this 
analysis: an electron (e) hypothesis, a muon (/i) hypoth- 
esis, a two-photon (77) hypothesis, and a pion (ir°) hy- 
pothesis. The electron and muon fits are single track 
fits parametrized by vertex position (x, y, z,i), direction 
(6,<j)), and energy (E). The probability of the charge 
and time of each PMT hit resulting from a given track 
configuration can be estimated using an optical model in- 
cluding predictions for Cherenkov and scintillation light 
emission profiles for the outgoing lepton and a descrip- 
tion of light propagation in the detector. The optical 
model is informed by in situ measurements. For each 
event, the negative-log-likelihood of the prediction com- 
pared to data is minimized over the space of track config- 
urations. The muon and electron hypotheses differ most 
significantly in the predicted topology of their associated 
Cherenkov rings. Rings from electrons are blurred by 
multiple scattering and electromagnetic showers whereas 
muons, with straighter tracks and no associated shower- 
ing, project sharp rings onto the surface of the detector. 
The two-photon hypothesis is a two-track fit. Concep- 
tually, the two tracks represent the two photons from a 
7T° decay. In practice, each track is treated using the 
electron hypothesis since photons resemble electrons in 
the detector. The two tracks share a common vertex and 
are parametrized by direction and energy as in the one- 
track fit and are each also parametrized by the photon 
conversion length. The ir° hypothesis is enforced by con- 
straining the photon-photon invariant mass m 77 to the 
7T° mass in the two-photon fit. Reconstructed variables 
are used to further refine the NC l7r° sample. We re- 
quire interaction vertices of candidates to (4) be within 
a 500 cm-radius fiducial volume according to the elec- 
tron fit. Candidates must favor the electron likelihood 
over the muon likelihood: more precisely, we require (5) 
\og{C e / Cfj) > 0.05. The distribution of this difference 
appears in Fig. [T] The separation between events with 
and without a tt is evident. Candidates must then fa- 
vor the pion likelihood over the electron likelihood: (6) 
\og(£ e /£ 7T ) < 0. Finally, we require that (7) the invari- 
ant mass extracted from the two-photon fit reside in the 
interval [80,200] MeV/c 2 . Figure Fy includes the invariant 
mass distribution; a distinct peak around the ir° mass 
of 134.97 MeV/c 2 is visible. Only a miniscule number of 
events in the mass peak is predicted to contain no 7r°s. A 
summary of the effect of each cut on the predicted purity 
and efficiency of each sample appears in Table [TTJ 

With 6.46 x 10 20 protons-on-target (POT) collected in 
neutrino mode running, 21375 events pass the selection 
requirements. In antineutrino mode running, 2789 events 
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FIG. 1: (a) Distribution of the difference between the e log- 
likelihood and the /i log- likelihood for events passing cuts (1)- 
(4) described in the text for neutrino mode running (top) 
and antineutrino mode running (bottom). Monte Carlo is de- 
picted by a dark-gray line and data by black dots. Both data 
and Monte Carlo are absolutely normalized to 10 20 POT. Er- 
ror bars are statistical only. Also shown are the contribu- 
tions from events containing no 7r° in the detector (translu- 
cent light-gray fill), signal NC lir° production (dark-gray fill), 
and incoherent (hatched fill) and coherent (gray fill) exclusive 
NC l7r° production according to identification at the neutrino 
interaction vertex. Candidate NC l7r° events are selected in 
the region indicated by the arrows, (b) Distribution of the 
reconstructed 7-7 invariant mass for events passing cuts (1)- 
(6) described in the text. The dashed vertical line marks the 
expected tt° mass. 



pass selection requirements with 3.68 x 10 20 POT col- 
lected. The Monte Carlo underestimates the number of 
events passing the cuts in neutrino mode by 10.9(8) stQt % 
and overestimates it in antineutrino mode by 5(2) stat %. 
In each running mode, the sample collected is the largest 
set of NC l7r° events recorded to date. These samples 
exceed the total of all samples collected by previous ex- 
periments by roughly an order of magnitude. 
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TABLE II: Predicted purity of the NC lvr sample and 
NC 1ty° selection efficiency in neutrino and antineutrino 
modes after each cut described in the text. Purity includ- 
ing wrong-sign induced signal sources is presented parenthet- 
ically. 



IV. ANALYSIS AND RESULTS 

A selection of photon kinematic distributions from 
the 7T° fit appears in Fig. [2J An incorrect prediction 
of 7r°s in the final state accounts for the disagreement 
between data and Monte Carlo in these distributions 
rather than any failure of the reconstruction, which has 
been separately vetted[3B]. Correcting the Monte Carlo 
with an in situ measurement of the rate of ir° produc- 
tion as a function of momentum — a kinematic that is 
strongly influenced by FSI — improves the level of agree- 
ment substantially [35] ■ The photon kinematics are used 
to derive the ir° kinematics. The four-momentum of the 
7r is simply the sum-momentum of the two photons. The 
incoming neutrino is assumed to be traveling in the beam 
direction, which is oriented with the z axis by convention, 
so the 7T° angle is taken to be the angle relative to the 
z axis. Using the partitions appearing in Fig. [3j we gen- 
erate histograms of ir° momentum and 7r° angle for the 
NC l7T° candidates. The neutrino mode 7r° momentum 
distribution extends to 1.5 GeV/c while the antineutrino 
mode distribution extends to 1.1 GeV/c. 

Background events arise from wrong-neutrino induced 
NC l7r° production and interactions in the detector mim- 
icking the signal signature. Interactions occurring out- 
side the detector ("dirt events") introduce negligible 
background. The fractional composition of the back- 
ground is listed in Table [TTTJ Of the wrong-neutrino back- 
grounds, only v^s in the beam constitute a significant 
background. Indeed, because of the sizable contamina- 
tion in the beam, wrong-sign production is the dominant 
background to the measurements; the measure- 
ments are relatively unaffected by wrong-sign production. 
The v e {ue) induced background is very small by virtue of 
the small beam contamination. The size of the detector 
affects the probability that particles emerging from the 
target nucleus will produce a n° in the tank. To avoid 
influencing the measurement with detector geometry, we 
include events with a tt° produced anywhere outside the 
target nucleus (and no ir° exiting the initial target nu- 
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FIG. 2: (a) The distribution of the reconstructed energy of the more energetic 7 from the 7r° decay in NC l7T° candidates from 
Monte Carlo (dark-gray line) and data (black dots). Results from neutrino mode running appear on the top and antineutrino 
mode running on the bottom. Error bars are statistical only and distributions are absolutely normalized to 10 20 POT. (b) 
The reconstructed energy of the less energetic 7. (c) The reconstructed opening angle between the two photons, (d) The 

reconstructed total momentum in the beam direction. 
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FIG. 3: (a) The reconstructed ir momentum distribution for 
NC l7T° candidates in neutrino mode running. The Monte 
Carlo distribution is shown as a dark-gray line and data as 
black dots. The box histogram is the systematic error on the 
Monte Carlo distribution; the error bars on the data are sta- 
tistical only. Distributions are absolutely normalized to 10 20 
POT. The black filled histogram is the non-NC 1tt° back- 
ground and the hatched histogram above is the additional 
contribution from wrong-sign induced NC 1tt° production, 
(b) The reconstructed 7r° angle distribution in neutrino mode 
running, (c) The reconstructed tv° momentum distribution in 
antineutrino mode running, (d) The reconstructed ty° angle 
distribution in antineutrino mode running. 



cleus) as background. Background interactions typically 
mimic signal events through a combination of the pro- 
duction of a 7T° outside the target nucleus and missed 
detection of other outgoing particles. NC tt ± production 
at the neutrino vertex is the most significant background 
to our signal. The w ± can readily charge exchange into 
a 7T°. NC elastic, multipion, CC ^ , and CC n° inter- 
actions each contribute to the background at a similar 
level. CC ir^~ events mimic the signal in the same man- 
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TABLE III: Predicted fractional composition of NC lw back- 
ground in neutrino and antineutrino modes broken down by 
exclusive channel at the initial neutrino interaction vertex and 
wrong-neutrino source. 



ner as their NC counterparts but also require that the 
outgoing lepton is undetected (captured or low momen- 
tum). NC elastic events contribute via ir° production 
induced in the detector by the outgoing nucleon and mul- 
tipion events through interactions producing a dominant 
7T°. FSI creating additional mesons cause a small frac- 
tion of incoherent exclusive NC l7r° events to be actually 
classified as background. 

As the initial step to extract the cross section, the 
Monte Carlo prediction of the background rates is used to 
extract the signal rate from the NC In sample. We sub- 
tract the absolutely normalized rate of all backgrounds 
except the wrong-sign NC l7r° background from the rate 
of candidate events in each bin of the kinematic distri- 
butions. To remove the wrong-sign content, we multiply 
the remaining content of each bin by the estimate of the 
right-sign NC l7r° fraction in that bin. 

Biases in the reconstruction, as well as detector ef- 
fects, smear the measured kinematics of the outgoing 
pion. This distortion is characterized in the response 
matrix, R. For a measurement, x, and a partition of 
the domain of x, (X n ), Rij is the probability that the 
reconstructed value of x is in bin i of (X n ) if the true 
value of x is in bin j. The response matrices for our four 
measurements, as estimated by Monte Carlo, appear in 
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FIG. 4: (a) Response matrix for the measurement of 7r° mo- 
mentum in NC 1tv° events satisfying selection cuts in neutrino 
mode, (b) The same for ty° angle in neutrino mode, (c) The 
same for ty° momentum in antineutrino mode, (d) The same 
for 7T° angle in antineutrino mode. 



Fig. [4] The response matrices indicate a tendency of the 
reconstruction to slightly overestimate ir° momentum, es- 
pecially at low momentum. In contrast, the response ma- 
trices for the measurement of tt° angle demonstrate little 
bias and excellent resolution in the forward region. In 
order to produce a physically meaningful measurement 
rather than one idiosyncratic to the experiment, we cor- 
rect the measurement for this distortion using a process 
known as unsmearing (or unfolding). Since the and 
Pfj, distributions differ in statistics by an order of mag- 
nitude and the p^o and cos v o distributions differ rad- 
ically in shape, using only one unsmearing technique is 
not necessarily appropriate. We evaluate three options — 
applying (1) Tikhonov regularized unsmearing with the 
regularization strength chosen by the SVD prescription 
detailed by Hocker and Kartvelishvili[37j, (2) a method 
analogous to one iteration of a Bayesian approach de- 
scribed by D'Agostini[35], and (3) no unsmearing — and 
select the least-biased result according to an unsmearing 
bias estimate from Cowan [39]. The unsmearing methods 
are described in greater detail in Appendix [5] We do not 
use matrix inversion to unsmear since it produces results 
with unacceptably large variance. We apply method (1) 
to the p^o distribution, method (2) to the cos^o 
and z/p p^o distributions, and method (3) to the cos 9^o 
distribution. 

After unsmearing the kinematic distributions, we ap- 
ply corrections to compensate for the misestimation of 
the number of events in the fiducial volume due to mis- 
reconstructed interaction vertices and losses due to de- 
tection inefficiency. These corrections appear in Fig. [5] 
In the former case, a bias in the reconstruction to pull 
interaction vertices to the center of the detector leads 
to a 7% excess of events being counted in the fiducial 
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FIG. 5: (a) NC In selection efficiency (solid line) and frac- 
tional decrease in the number of events in the fiducial volume 
when using the true vertex versus the reconstructed vertex 
(dashed line) as functions of n° momentum in neutrino mode, 
(b) The same for tv° angle in neutrino mode, (c) The same 
for 7T° momentum in antineutrino mode, (d) The same for ty° 
angle in antineutrino mode. 



volume. We subtract the fraction of nonfiducial events 
from each bin of each distribution. The average NC l7r° 
selection efficiency for each measurement is 36%. The 
selection efficiency is momentum dependent: it is dimin- 
ished at high and low momentum. At low momentum, 
the log(£ e /£ /J ) cut becomes more inefficient as the ability 
of the reconstruction to discriminate between muonlike 
and electronlike events is reduced. At higher momentum, 
loss of containment causes a larger proportion of signal 
events to fail the veto PMT hits requirement. Loss of con- 
tainment is responsible for the rejection of 11% of signal 
events in neutrino mode and 13% in antineutrino mode. 
To recover the rate of events, we divide the kinematic 
distributions by the efficiency in each bin. 

With the rate of NC l7r° production recovered, we 
must divide by the integrated flux and the number of 
targets to recover the flux-averaged cross section. We 
predict the flux at MiniBooNE using a GEANT4-based 
simulation of the neutrino beam [31]. Primary inter- 
actions of beam protons on the Be target producing 
7r ± s, K 0,± s, protons, or neutrons are handled by a cus- 
tomized framework incorporating external data. In par- 
ticular, the prediction of charged pion production (which 
is the dominant source of f„ and z/„) is based on data 
from HARP gD] and BNL E910gT]. The flux predic- 
tion in both neutrino and antineutrino modes appears 
in Fig. [6] The simulation predicts an integrated flux of 
(3.35 ± 0A3 sys ) x 10 11 v^/cm 2 over the course of neu- 
trino mode running and (1.08 ± 0.12 sys ) x 10 11 v^/cm 2 
over antineutrino mode running. The uncertainty in the 
flux in neutrino (antineutrino) running can be split into 
12.1% (13.1%) from secondary meson production uncer- 
tainties, 4.1% (2.8%) from the horn magnetic field (skin 
depth and current variations) and secondary interactions 
outside of the target, and 2% (2%) from the accounting 
of the number of protons delivered on target. Using a 
measured value of 0.845 ± 0.001 g/cm 2 for the density of 
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FIG. 6: The predicted flux of z/ M (solid lines) and z5 M (dotted 
lines) at the MiniBooNE detector with the horn configured 
in neutrino mode (black lines) and antineutrino mode (gray 
lines). The flux prediction is available at the MiniBooNE 
website HU. 



the mineral oil in the detector, we can determine that 
there are 2.664 ± 0.003 x 10 32 nucleons in the 500 cm- 
radius fiducial volume. Dividing each differential rate by 
the number of targets and the appropriate integrated flux 
yields the flux-averaged cross section per nucleoli. 

Plots of the resulting absolute differential cross sec- 
tions for NC l7r° production on CH2 appear in Fig. [7] 
and the tables in Appendix [C] Per our signal defini- 
tion, these cross sections include the effects of final state 
interactions. Integrating the differential cross sections 
yields total cross sections of (4.76 ± 0.05 stat ± 0.76 sys ) x 
10~ 40 cm 2 /nucleon at a mean energy of (E u ) = 808 MeV 
for f M -induced production and (1.48±0.05 s t a t±0.23 sl , s ) x 
10~ 40 cm 2 /nucleon at a mean energy of (E u ) = 664 MeV 
for i/^-induced production. These cross sections are flux- 
averaged; hence, they are specific to the neutrino flux at 
MiniBooNE [42 . Being the first absolute measurements 
of NC l7r° production, there are no other measurements 
with which to compare. 



SYSTEMATIC UNCERTAINTIES 



gle pion, coherent single pion, and multipion production 
are varied by 6.2%, 25%, 27%, and 40% from their cen- 
tral values of 1.23 GeV/c 2 , 1.10 GeV/c 2 , 1.03 GeV/c 2 , 
and 1.30 GeV/c 2 , respectively. The binding energy and 
Fermi momentum values used in the relativistic Fermi gas 
model[44 underlying the simulation of QE, NC elastic, 
and incoherent pion production are varied by 26% and 
14% from their central values of 34 MeV and 220 MeV/c, 
respectively. The total normalization of QE scattering, 
deep inelastic scattering (DIS), and A radiative processes 
are varied by 10%, 25%, and 12.2%, respectively. A Pauli 
blocking scale factor for CC QE events, k 45J, is varied 
by 0.022 from its central value of 1.022. In the target nu- 
cleus, the cross sections for pion absorption, pion charge 
exchange, and A interactions (AN — > N'N), are varied 
by 25%, 30%, and 100%, respectively. Pion scattering 
cross sections in the mineral oil outside the target nu- 
cleus are varied by 35% for absorption and 50% for charge 
exchange. The uncertainty in our pion interaction simu- 
lation is validated using external data for interactions on 
carbon [46Tl49] , In total, cross section uncertainties con- 
tribute an 8.4% uncertainty in the measured v„ NC l7r° 
production cross sections and 7.7% in the cross sec- 
tions. 

Uncertainty in the optical model in the detector and 
PMT response as well as bias in the unsmearing make up 
the detector uncertainties. Optical model uncertainties 
include variations in the amount of light production and 
in the propagation of light in the detector. A total of 39 
parameters are varied. For the PMT response, we assess 
one uncertainty by adjusting the discriminator threshold 
in the DAQ simulation from 0.1 PE to 0.2 PE and another 
by generating an excursion in the charge-time correlation 
of PMT hits. We also assess the estimated bias in the 
unsmearing as an error. Since unsmearing preserves the 
number of events in a distribution by design, the bias 
produces only a small uncertainty on the normalization 
of the cross section; the error is principally in the shape. 
Detector uncertainties constitute a 5.1% uncertainty in 



the cross section and 4. 



% in the cross section. 



Systematic uncertainties can be grouped into three 
principal categories — flux related, cross section related, 
and detector related. We gauge the uncertainty in the 
measurements including bin-to-bin correlations by calcu- 
lating the covariance of the measurements over a set of 
Monte Carlo excursions wherein underlying parameters 
are varied within their uncertainties and correlations. 

The same uncertainties affecting the integrated flux 



prediction detailed in Sec. IV also affect the Monte Carlo 
predictions used in the cross section calculation, e.g. the 
background prediction. In total, flux uncertainties pro- 
duce a 12.4% overall uncertainty in the cross sections 
and 12.7% in the cross sections. 

The cross sections associated with background pro- 
cesses are varied within their uncertainties. The rele- 
vant axial masses for quasielastic (QE), incoherent sin- 



VI. DISCUSSION 

Honing models of single pion production continues 
to be of theoretical interest. In particular, elucidating 
the nature of coherent pion production is a very active 
pursuit [31 141 [5H51 [TUrUT] . As an illustration, our own pre- 
diction of single 7r° production can be tested against our 
data. 

We predict single 7r° production using models by Rein 
and Sehgal[2j [5] as implemented in Nuance. The axial 
masses for incoherent and coherent pion production are 
assumed to be 1.1 GeV/c 2 and 1.03 GeV/c 2 , respectively. 
Additionally, we use the Nuance FSI simulation in lieu 
of the pion absorptive factor suggested by R-S for co- 
herent pion production. Assuming these predictions [5 1). 
MiniBooNE found that coherent pion production com- 
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FIG. 7: Flux- averaged absolute differential cross sections for NC 1tv° production on CH2 including the effects of FSI. Data are 
shown as black dots with statistical error bars and systematic error boxes. The dark-gray line is the Monte Carlo prediction [33] 



using R-S models of single pion production 2, 5 modified as described in the text, (a) 



dp 



for ^-induced production, (b) 



for ^-induced production, (c) 



dp 



for ^-induced production, (d) 



for f M -induced production. The numerical 



values for the cross sections appear in Appendix [c] and are also available at the MiniBooNE website [43] . 
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FIG. 8: (a) The differential cross section for inclusive NC i7r° 
production as a function of tv° angle in neutrino mode (above) 
and antineutrino mode (below). Data is indicated by black 
dots with statistical error bars and systematic error boxes. 
The Monte Carlo prediction including the R-S single pion 
production models [2l [5] as implemented in Nuance adjusted 
according to [35] is indicated by the thick black line. The 
prediction omitting the coherent portion of NC 7r° production 
is indicated by the dashed black line. The arrow indicates the 
region for which a \ 2 is quoted in the text. The horizontal 
scale is magnified in the forward region. 



prises (19.5 ± l.l sta t ± 2.5 sys )% of exclusive NC lir° pro- 
duction in neutrino mode [35 ]. This fraction implies a 
35% reduction in R-S coherent pion production (and a 
corresponding 5% increase in incoherent production) that 
is incorporated into our Monte Carlo prediction. Figure[8] 
compares the differential cross section in tt° angle (the 
distribution most sensitive to the production mode) from 
data to our Monte Carlo prediction with and without 
coherent pion production. In the forward region above 
cos 9^0 = 0.6, the x 2 between neutrino (antineutrino) 
data and the Monte Carlo including coherent pion pro- 
duction is 8.23 (13.6) with 9 (5) degrees of freedom, which 
corresponds to a p- value of 0.511 (0.018). Without coher- 
ent pion production, the \ 2 worsens to 45.1 (25.7) with 9 
(5) degrees of freedom, which corresponds to a p- value of 
8.7 x 10~ 7 (0.0001). Both the neutrino and antineutrino 
data clearly favor the model of single ir° production with 
nonzero coherent content. Though the model including 
coherent pion production is favored, the shape disagree- 
ment evident in Fig. [8] substantiates, but does not con- 
firm, the claims[HG>] that the R-S model[2] is inadequate 
at neutrino energies below 2 GeV. Alternative mecha- 
nisms, such as an incorrect prediction of the FSI[52], can 
account for the disagreement in part, but they are un- 
likely to explain the discrepancy in full, particularly in 
antineutrino mode. Used in concert, our measurements 
in momentum and angle can be used to evaluate and re- 
fine the abundance of modern models that endeavor to 
correctly describe single pion production on nuclei with 
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the effects of other mechanisms disentangled. 



Our measurement is designed to be independent of the 
assumed models of single pion production and FSI. Al- 
though, in making a pure or measurement with 
a contaminated beam, we introduce some dependence on 
the assumed single pion production model by subtracting 
wrong-sign content. In Appendix [XJ we characterize this 
sensitivity and present an alternative, fully-independent 
measurement. 



In addition, we assess the cross section for and P„ 
induced incoherent NC l7r° production defined at the ini- 
tial neutrino interaction vertex as a means to compare 
with past measurements. Such an exclusive measure- 
ment is naturally quite sensitive to assumed models of 
both single pion production and FSI. We use the same 
selection cuts as in the primary analysis. Because co- 
herent NC lir° production is a background to this mea- 
surement, the result suffers from a fairly low predicted 
signal fraction: 57% in neutrino mode and 34% in an- 
tineutrino mode. We use the same selection of unsmear- 
ing techniques used in the primary analysis as well. The 
nonfiducial fraction is also predicted to be the same at 
7%. Unlike in the inclusive measurement, the efficiency 
correction includes a correction for FSI predicted using 
Monte Carlo that recovers the kinematic distributions 
at the initial neutrino interaction vertex. This over- 
all efficiency including selection inefficiency and FSI is 
predicted to be 24% in both neutrino and antineutrino 
modes. After all corrections, we find the cross section to 
be (5.71±0.08 stat ±1.45 sys ) x 1CT 40 cm 2 /nucleon for v„- 
induced incoherent exclusive NC l7r° production on CH2 
and (1.28 ± 0.07 stat ± 0.35 sys ) x 1CT 40 cm 2 /micleon for 
i^-induced production. These cross sections are averaged 
over the MiniBooNE flux as well. Here, the significance 
of FSI becomes apparent: the incoherent exclusive 
NC l7r° production cross section actually exceeds the 
inclusive NC l7r° production cross section. Repeating the 
measurement using the models of 0] and [3] discussed in 
Appendix [A| yields values of (6.51 ± 0.08 stat ± 1.56 SJ/S ) x 
1(T 40 cmTnucleon and (6.20 ± 0.08 stat ± 1.52 sys ) x 
10~ 4n cm 2 /nucleon, respectively, for induced produc- 
tion, and (1.78±0.07 stat ±0.42 sys ) x 1CT 40 cm 2 /nucleon 
and (1.62 ± 0.07 sfat ± 0.39 aJ , a ) x 1CT 40 cm 2 /nucleon, re- 
spectively, for D„ induced production. The variation 
in the measurements extracted under alternative models 
of coherent pion production illustrate the model depen- 
dence of the extracted incoherent cross section. These 
measurements are plotted against prior measurements 
and the Nuance prediction (using R-S) in Fig. [9] A 
comparison can be made only to the result of the reanal- 
ysis of the Gargamelle data[24] since the measurement at 
Aachen-Padova was limited to production on protons |23j. 
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FIG. 9: (a) The flux-averaged total cross sections for v^- 
induced incoherent exclusive NC lw production on CH2 cor- 
rected for FSI. Points 1, 2, and 3, are the cross sections 
extracted using the MiniBooNE implementation of the R-S 
model for coherent pion production, the model in 0], and the 
model in [3], respectively. The points are placed at the mean 
energy of the beam in neutrino mode; the spread is only for 
clarity. The curve is the Nuance prediction using the R-S 
model. Also shown for comparison is the measurement made 
from the Gargamelle data 24 . The Gargamelle experiment 
used a propane and freon (C3H8 + CFaBr) target, (b) The 
same for z/ M -induced incoherent exclusive NC l7r° production. 
In this case, there are no external measurements to compare 

to. 



VII. CONCLUSION 

In conclusion, we have used the largest sample of NC 
l7T° events collected to date to produce measurements of 
absolute differential cross sections of NC l7r° production 
induced by both neutrinos and antincutrinos on CH2 as 
functions of both ir° momentum and ir° angle averaged 
over the MiniBooNE flux. These measurements, which 
are the principal result of this work, can be found in Fig. [7] 
and Table IIV1 The total cross sections have been mea- 
sured to be (4.76±0.05 stot ±0.76 sas ) x 1CT 40 cm 2 /nucleon 
for interactions at a mean energy of 808 MeV and 
(1.47 ± 0.05 stQt ± 0.23 sys ) x 10~ 40 cm 2 /nucleon for 
interactions at a mean energy of 664 MeV. These mea- 
surements should prove useful to both future oscillation 
experiments seeking to constrain their backgrounds and 
those developing models of single pion production seek- 
ing to test their predictions. We have additionally mea- 
sured total cross sections for incoherent exclusive NC l7r° 
production on CH2 to compare to a prior measurement. 
These cross sections were found to be (5.71 ± 0.08 sta t ± 
1.45 sys ) x 10 -40 cm 2 /nucleon for ^-induced production 
and (1.28 ± 0.07 stat ± 0.35 sys ) x 10~ 40 cm 2 /nucleon v^- 
induced production. 
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FIG. 10: Pion angular distributions in the forward region for 
NC l7r° candidates in (a) neutrino and (b) antineutrino mode 
running. Data is indicated by black dots with statistical er- 
ror bars. The Monte Carlo prediction using the reseated [35] 
R-S model of coherent pion production as implemented in 
Nuance[21 [33] is indicated by the solid black line with gray 
systematic error boxes. The predictions using the models of 
H] and [3] are indicated by the dotted line and the dashed 
line, respectively. The systematic error in the predictions us- 
ing the alternative models is of the same relative size as the 
prediction using R-S; it is omitted for clarity. Distributions 



are normalized to 10 



POT (b) 
mode. 



The same for antineutrino 



Appendix A: Measurement Model Dependence 

Subtraction of wrong-sign induced NC l7r° signal 
events inevitably couples our measurements to the as- 
sumed model of NC lir° production. For the sake of 
example, we considered the effect of substituting the co- 
herent pion production models of Refs. [4] & [3J into our 
Monte Carlo prediction. The difference in the angular 
distribution of events satisfying the NC l7r° selection cuts 
under these models appears in Fig. 10 Both the micro- 



scopic models demonstrate a sharper peaking in forward 
direction compared to the MiniBooNE R-S central value. 
However, owing to a different choice for the N-A transi- 
tion axial form factor C^, Ref. [4] predicts substantially 



less production than Ref. [3J. In Fig. 11 the ratio of the 
angular cross sections extracted assuming the models in 
Refs. [31 0] relative to the primary result is shown. Be- 
cause of the low wrong-sign contamination, the cross 
section is relatively insensitive to changes in the model; 
however the cross section deviates more significantly 
under the model variations. The total cross section de- 
creases by 5.8% under [4] and 4.4% under [3j; the total 
cross section varies by < 1% in either case. Even though 
an attempt is made to partially mitigate model depen- 
dence in the wrong-sign subtraction by scaling by the 
right-sign fraction rather than outright subtracting the 
rate, the large wrong-sign fraction in antineutrino mode 
together with the very large variation from [4 conspire to 
generate a non-negligible difference in the measured cross 
section. Such dependence is unavoidable when measuring 
a D^-only cross section. 

In order to provide a measurement that is unbiased 
by any assumed model of NC l7r° production, against 
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FIG. 11: (a) The ratio of the NC lw angular cross section 
extracted using the coherent production models in [3] (dotted 
line) and [3] (dashed line) compared to the principal result 
extracted using the R-S model, (b) The same for the 
NC l7r° cross section. 




FIG. 12: Flux-averaged absolute differential cross sections for 
fM+^M -mcra ced NC lit production on CH2. Data are shown 
as black dots with statistical error bars and systematic error 
boxes. The dark-gray line is the Monte Carlo[33] prediction 
using R-S models of single pion production]!] [5] modified as 
described in the text, (a) , in neutrino mode, (b) , a — 

in neutrino mode, (c) , drT in antineutrino mode, (d) , d<T a — 

V ' dp q V / a cos u 

in antineutrino mode. 



which other models can be tested, we performed the 
principal analysis again in the exact same manner ex- 
cept signal events induced by wrong-sign neutrinos are 
not subtracted. These combined v^+v^ measurements 
are almost entirely free of the model dependence intro- 
duced by the wrong-sign subtraction at the cost of be- 
ing a less immediately meaningful measurement. Nat- 
urally, the signal fraction increases: it is 75% in neu- 
trino mode and 82% in antineutrino mode. The non- 
fiducial fraction and selection efficiency remain the same 
(7% and 36%, respectively). The combined integrated 
flux over neutrino mode running is (3.57 ± 0.50 sys ) x 
10 11 (f fl +^)/cm 2 and the combined integrated flux over 
antineutrino mode running is (\.'b?>±0.2\ S y S ) x 10 11 + 
C" M )/cm 2 . We find the flux-averaged total cross section 
for -(-^-induced NC l7r° production on CH 2 to be 
(4.56±0.05 stat ±0.71 sys ) x 10~ 40 cm 2 /nucleo n in neutrino 
mode and (1.75±0.04 sto t ±0.24 sys ) x 10~ 40 cm 2 /nucleon 
in antineutrino mode. The v^+v^ differential cross sec- 
tions appear in Fig. [12] 
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Appendix B: Unsmearing 

Wc begin by defining an abstract unsmearing scenario. 
Suppose we make a measurement of a variable x over an 
n-bin partition of the domain of x, (X n ), that is sub- 
ject to smearing dictated by a response matrix R. If 
the discrete probability density function for x over the 
partition is a, then the probability distribution function 
for the measured values is (3 = Ra. In an actual mea- 
surement of N events, we make a draw b ~ N/3 which 
corresponds to an unknown true distribution a ~ Na. 
In unsmearing, we seek to determine an estimator for a, 
a, knowing only b and Monte Carlo estimates of R and 
a, R MC and ot MC . In this analysis, we treat smearing as 
affecting only the shape of a distribution and not the nor- 
malization as including efficiency losses would do. Here 
we describe three unsmearing methods, two of which are 
used in the analysis. 

A naive method of unsmearing follows from the expres- 
sion (3 = Ra given the population distributions and the 
response matrix. It follows that a. = R _1 /3. Hence, if 
R MC estimates R well, then we may choose 



a = R 



MC~ 



(Bl) 



to be an estimator for a. This choice of unsmearing is 
knows as matrix inversion. Since Eq. |B1| involves the 
inversion of a matrix, it is particularly sensitive to per- 
turbations in R MC and b. Matrix inversion often proves 
to be too unstable to be useful. 

The second method is a specialization of Tikhonov reg- 
ularization. Under Tikhonov regularization we choose 
the a that minimizes the quantity 



(R A/C a - b) T V (b) (R MC a - b) + r||La| 



(B2) 



where V (b) is the covariance matrix for b, L is some lin- 
ear operator, and r is a constant controlling the strength 
of regularization. The quantity on the left is simply a \ 2 
between the measured reconstructed distribution and the 
smeared estimator for the true distribution. Minimizing 
only the \ 2 results in the estimator a = R MC 1 b — the 
result of matrix inversion. This result is usually highly 
unstable. The right hand term is a regularizing term that 
reduces the variance by adding a penalty for not satis- 
fying some a priori characteristic of a encoded by the 



action of L. For this analysis, we assume that the true 
distributions are smooth, so we seek to minimize the cur- 
vature of the estimate. To that end, we choose L to be the 
second finite-difference operator (a discretization of the 
second derivative). Equation B2 can by minimized ana- 
lytically. Typically no constraint is placed on the mini- 
mization, but we use the method of Lagrange multipliers 
to minimize under the constraint that a, = J^. bi per 
our objective to not change the normalization, which re- 
sults in 




R MC + rVR Mcr L T L 



J2jkUjiVl m R M kr, 



(B3) 



The choice of r follows the prescription in Ref. |37j . 
Bias is introduced not through the Monte Carlo, but the 
choice of Tikhonov matrix, L. 

The third method is equivalent to a single iteration 
of the Bayesian method described in Ref. [35]. Since 



MC 



3 J % 



1 V i by definition, it follows that S 



MC 1 



(1, 1, 1) = (1, 1, 1). We construct a matrix U, given 
by 



U = diag(a MC )S MC ' diag(/3 iuc; ) 



MC\-1 



(B4) 



By construction U/3 



MC 



Assuming that the 



Monte Carlo is a good estimator for the data, then we 
can use a = Ub as an estimator for a. This method 
introduces bias from the Monte Carlo. 



Appendix C: Cross Section Values 



The Vn- and ^-induced NC l7r° production cross sec- 



tion measurements on CH2 are tabulated in Table |IV| 
The measurements together with full error matrices 
are also in a data release available at the MiniBooNE 
website H3] . 
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Plr o (GcV/c) (io 


da 1 dp o 

- 39 cm 2 /(GcV/c)) 


Pw o (GcV/c) (iQ 


da /dp q 

- 39 cm 2 /(GeV/c» 


(0.00,0.10) 
(0.10,0.15) 
(0.15,0.20) 
(0.20,0.25) 
(0.25,0.30) 
(0.30,0.40) 


0.18 ± 0.06 
1.19 ± 0.21 
1.63 ± 0.24 
1.58 ± 0.21 
1.28 ± 0.21 
0.87 ± 0.14 


(0.40,0.50) 
(0.50,0.60) 
(0.60,0.80) 
(0.80, 1.00) 
(1.00, 1.50) 


0.47+ 0.09 
0.21 ± 0.06 
0.05 + 0.04 
0.03+ 0.02 
0.01 ± 0.01 


(a) Vfj, NC 1tt° production p^o differential cross section 


cos 9 o 


da jd cos 9 o 
(10- 40 cm 2 /l) 


cos 9 o 


da 1 d cos 9 o 
(10- 40 cm 2 /l) 


(-1.000, -0.620) 
(-0.620, -0.340) 
(-0.340, -0.130) 
(-0.130, +0.060) 
(+0.060, +0.200) 
(+0.200, +0.320) 
(+0.320, +0.420) 
(+0.420, +0.520) 
(+0.520, +0.600) 


0.82 ± 0.16 
1.10 ± 0.22 
1.37 + 0.24 
1.64+ 0.27 
1.99 ± 0.43 
2.26 ± 0.35 
2.58 ± 0.43 
2.82 ± 0.46 
3.16 ± 0.50 


(+0.600, +0.670) 
(+0.670, +0.730) 
(+0.730, +0.780) 
(+0.780, +0.830) 
(+0.830, +0.870) 
(+0.870, +0.910) 
(+0.910, +0.950) 
(+0.950, +0.975) 
(+0.975, +1.000) 


3.68 + 0.63 
3.94 + 0.65 
4.38+ 0.70 
4.96+ 0.77 
5.49 ± 0.91 
6.33+ 1.04 
7.28+ 1.18 
8.42 ± 1.45 
9.56 ± 1.58 


(b) NC lit production cos O^o 


differential cross section 


P x o (GcV/c) <io 


da/dp^o 

- 40 cm 2 /(GoV/c)) 


P w o (GcV/c) (iq 


da/dp w o 

- 40 om 2 /(GeV/c)) 


(0.00,0.13) 
(0.13,0.17) 
(0.17,0.21) 
(0.21, 0.24) 
(0.24,0.28) 


1.13 ± 0.25 
5.20 + 0.86 
5.86 ± 0.86 
5.26 ± 0.78 
4.42 ± 0.64 


(0.28,0.32) 
(0.32,0.37) 
(0.37,0.44) 
(0.44,0.57) 
(0.57, 1.10) 


3.68 + 0.55 
2.84 + 0.49 
1.72 ± 0.36 
0.71 ± 0.19 
0.11 ± 0.06 


(c) Uft NC l7r° production p^o differential cross section 


cos 6^o 


da / d cos 9^0 
(10- 40 cm 2 /l) 


cos 9^o 


da j d cos 9 o 
(10- 40 cm 2 /l) 


(-1.00, -0.60) 
(-0.60, -0.22) 
(-0.22, +0.12) 
(+0.12, +0.40) 
(+0.40, +0.60) 


0.38 ± 0.08 
0.40 ± 0.08 
0.50 ± 0.10 
0.61 ± 0.12 
0.69 ± 0.15 


(+0.60, +0.74) 
(+0.74, +0.85) 
(+0.85, +0.91) 
(+0.91, +0.96) 
(+0.96, +1.00) 


1.00 ± 0.20 
1.33 + 0.27 
1.94 + 0.38 
2.76 + 0.50 
4.06 ± 0.74 



(d) Vfj, NC 1tt° production cos 9^o differential cross section 

TABLE IV: Tabulated values of the flux-averaged differential 
cross sections for f M - and ^-induced NC 1ty° production on 
CH2 corresponding to the plots in Figure [7] The error quoted 
with the cross section values is the quadrature sum of the 
diagonal statistical and systematic error. 
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